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Abstract 

Recent observations in the OPERA experiment suggest that the neutrino could 
propagate with speed that is supcrluminal. Based on early theoretical work on tachy- 
onic fermions we shall study a modification of the Gross-Neveu model in two dimen- 
sions. We shall see that the theory results to the dynamical generation of real and 
imaginary masses. These imaginary masses indicate the possibility that tachyonic 
solutions (or instabilities) could exist in the theory. The implications of a tachy- 
onic neutrino coming from astrophysical sources are critically discussed. Moreover, 
we present a toy model that consists of an {7(2, 2) invariant Dirac Lagrangian. This 
theory can have tachyonic masses as solutions. A natural mass splitting between the 
solutions is a natural outcome of the formalism. 
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Modified 2d Gross-Neveu Model and Extrema of the Effec- 
tive Potential 



The Gross-Neveu model describes Dirac fermions in 1+1 dimensions, that have four 
fermions interactions, thus being a very elegant and useful tool for studies of the strong 
interactions. It is based on non-perturbative techniques, namely the large N expansion 
and it results to the dynamical generation of a fermionic mass, a phenomenon known as 
dimensional transmutation. It is described by the Lagrangian density 

N N 2 

£ = J2 Wi-fWi + ( E Mi) (1) 

i i 

for i = 1, 2, ...N. The Lagrangian (JTJ is invariant under the chiral symmetry ipi —¥ 75^1- 
This symmetry is a kind of Z2 symmetry, since 7! = 1. Due to one- loop quantum correc- 
tions, this symmetry is broken by the theory. 

Recently an experimental result has put into question the nature of the neutrino. 
Particularly, the OPERA experiment [T] resulted that the neutrino propagates in space 
with speed that exceeds the speed of light. It is obvious that the experimental results 
must be re-examined in order to be sure that this result is valid. Nevertheless it is an 
intriguing result that deserves some attention and if proven true, much of the theoretical 
work involving the neutrino must be put in a new conceptual basis. Actually it has draw 
the attention of many theorists [2Hl6j . 

Insightful studies on the tachyonic nature of the neutrino was provided in the 80's by A. 
Chodos et. al [T5H20] . The authors proposed experiments that could prove the tachyonic 
nature of the neutrino field. They actually proposed the time of flight experiment, which 
is an OPERA-like experiment. Up to now, several papers followed, studying the tachyonic 
neutrino [21H30] . The Lagrangian that describes the tachyonic neutrino is given by [18) : 

CJ = ^757^(9^ — rmpil) (2) 

In the massless case, at the Standard Model level, the terms ^757^^^ an d iV>7 M c^V> are 
indistinguishable. 

In this brief letter we shall modify the Gross-Neveu Lagrangian according to the mass- 
less limit of the tachyonic Lagrangian ([2D - The modified Gross-Neveu Lagrangian is: 

N I N 2 

i i 

We shall dwell on the mass generation issues, and how the results are modified when using 
the above Lagrangian. We must note that the massless field theory of the tachyons does 
not suffer from the conceptual problems that the massive case has. Indeed no consistent 
Quantum field theory of tachyons exists up to now [IS], however it could be a good 
starting point to work towards the consistent incorporation of tachyons in the successful 
field theories. The Lagrangian ([3D is invariant under the chiral transformation tpi — > 75^. 
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We now compute the effective potential corresponding to ([3|). Using very well known 
techniques [31] we have: 
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where a is a composite field in terms of which we shall express the final expression of the 
effective potential. In the above, \x is a free parameter to be used for dimensional reasons 
and D = 2 — 2e. After performing the Gaussian integration over the fermions and taking 
the e — > we obtain: 
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Using the MS renormalization scheme to subtract the poles, the effective potential reads: 
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Minimizing the effective potential in respect to the parameter a, we obtain the equation: 

= (7) 
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The above equation can have four real roots of the form (a?i,— x\) and (2:2,— X2), with 
xi,X2 real positive number. This result is similar to the usual Gross-Neveu result. This 
means that a can take four real values, a = ±xi and a = ±X2- The field a is the minimum 
of the effective potential, thus the quantum theory can have two equivalent ground states. 
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Figure 1: The effective potential for g 2 N/ \i ~ 10. In the left figure we can see the minima 
of the potential at o\ = 100 and o<i = —100. In the right figure we can see the maxima of 
the potential 0-3 = 0.0091, 0-4 = -0.0091 

The fermions acquire finite masses which can be m = \a\ = x\ and m = \a\ = X2- It 
is obvious that the initial Z2 chiral symmetry is broken. Additionally, the fermionic fields 
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can have two different values for their mass, which is dynamically generated. The modified 
Gross-Neveu model we studied results to two equal masses (minima of the potential), a 
result similar to the original Gross-Neveu result. But the most intriguing difference in 
comparison to the Gross-Neveu is the existence of complex extremal points of the effective 
potential. Indeed a numerical study of the above equation ([7]) shows that there are values 
of the expression for which the equation has two complex roots. These roots indicate 
instability of the quantum theory, since the mass-square of the field a can be negative. 
But in view of the new experimental data we could argue that the above model Q can 
generate dynamically tachyonic masses. Recalling that the modified version of the Gross- 
Neveu we studied is based on a tachyonic Lagrangian, the fact that solutions for tachyonic 
masses exist, can be very useful. 

Let us give here some quantitative details of the calculation we performed. Recall that 
the original two dimensional Gross-Neveu model was described by the effective potential 

The minimization process results to the following equation: 

1 N a 2 

—r + —hi- 

2g 2 Air fj, 



2 ■ A „ - „2 " ( 9 ) 



It is obvious that equation ([9]) has only two real roots. In addition, if one tries to search 
for roots of the form a = ix, this results to mathematical inconsistencies. Hence no 
imaginary roots exists for the original two dimensional Gross-Neveu model. In the case 
of d?]), we have a plethora of roots. The most interesting are these of the form a = ix. 
Indeed, one can easily check numerically that if we put a = ix, then we wave two real 
solutions for x, namely x and —x. Hence the result is that, equation ([7]) has two purely 
imaginary roots, a = ix and a = —ix For example, if Ng 2 = 10 12 (large N-strong coupling 
limit) and fi = 10 6 , then o\ = 10 _7 i and 02 = —10~ 7 i. Therefore we have two equivalent 
tachyonic mass vacuum state solutions. Apart from purely imaginary solutions, equation 
([7]) has complex roots of the form a = a + bi, and as we mentioned earlier real roots. Let 
us discuss the implications of the complex roots at this point. Complex masses are very 
closely related to tachyonic fermion solutions. Indeed, a generalized form of Lagrangian 
©, found in [59] , results in the modified Dirac equation 



(i-fdp - A T )V = (10) 

with At = ia + 675. Under certain conditions, depending on the a and b choices, one can 
obtain tachyonic solutions in this case, but there are two more solutions. The authors 
of [29j call these "luxons" and "bradyons". To us, these solutions are under serious inves- 
tigation, since the solutions lead to serious quantization problems, more severe than the 
tachyonic solutions. In any case, if superluminal propagation in spacetime is true, then 
we should work and find a consistent quantization procedure, even by generalizing the 
theoretical formalism of quantization. As for the real roots that equation Q has, the first 
set zLx\ are local maxima and the other set ±2:2 axe local minima. Hence the situation is 
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equivalent to the 2d Gross-Neveu model, except the two maximum solutions. In figures 
1 (a) (left) and (b) (right) we plot the effective potential (|6|), for g 2 N/fi ~ 10. We can 
clearly see the two minima, o\ = 100, 02 = —100 and the two maxima 03 = 0.0091, 
04 = —0.0091. In figure 2 we have made a plot of the roots corresponding to equation (|Tj). 
for various values of g 2 N j \i. The y-axis represents the imaginary axis, while the x-axis, 
the real axis. As we can see there are purely imaginary, complex and real roots. 
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Figure 2: The minima of the effective potential V e ff f° r various values of the parameters 
g, ji in the large N limit 



Before closing this investigation, let us study the effective potential in various limiting 
cases of the parameters. We start with the case t < 1. Upon using the approximation 
for the logarithm, valid for x > 0: 



lnx = 2V r^-A—f" (11) 
^ (2k- 1) \x + lJ v ' 

and keeping the first term in the sum, the effective potential V e ff of equation © can be 
cast as: 

1 2 N A /a 4 
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If we minimize equation fjl2[) in respect to a , we obtain the following solutions: 

1 / q 2 Nm 2 A \ 
^ = ^+^ + 31!?. ( 13 ) 



3ov a 5 2 iv;' 
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60 V A g 2 N 

- s = i)( 2ra+i — f — + f N 

where we introduced the parameter A, which is equal to: 



A = [ m 3 N 3 g 6 - 2700mN 2 irg A + 30Vfa\ g 8 m 2 N 4 ( 1350vr - g 2 m 2 N) (14) 
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In the case we have ^ > 1, we cannot have analytic solutions. Nevertheless by using the 

fj, 

approximation, 

lnx = -K (x) -7 + ln2 (15) 

and in addition, the asymptotic form of the modified Bessel function Ko(x), valid for large 
values of x: 

Ko(*) ~ J^e-* (16) 

we obtain a simpler expression for the effective potential. Indeed the potential can be 

written as: 

1 2 N ,/ , fn ' 

The original Gross-Neveu model was renormalization scale invariant. We must say 
that the effective potential of the modified Gross-Neveu is also //-independent, as can be 
easily seen from the following equation 



v efA°) = T2 a2 + — a2 \^\h e ^-7 + ln2-2) (17) 



We shall not pursuit further the renormalization group issues in this paper. 



A Discussion on Tachyonic Fermions and Their Impact on Physical Phenomena 

In this section we discuss some issues in reference to the possible tachyonic nature of 
the neutrino. Obviously, the experimental results must be further scrutinized in order 
to be sure for their validity. However if true, the tachyonic nature of the neutrino could 
alter some aspects of physical processes or even give answers to unexplained mysteries. 
For example in hot-dense backgrounds, the plasmon decay and the effective neutrino- 
photon coupling must be reexamined in view of the tachyonic nature of the neutrino. 
The modification of the effective photon-neutrino interactions at high temperature could 
modify the neutrino emission from stars, or the neutrino emission from Gamma Ray 



6 



Bursts. On that account, an interesting impact of superluminal neutrino's is studied 
in [2]. In regards to Gamma Ray bursts, phenomenological results of quantum gravity 
predict that neutrino's could arrive faster than the photons |32j . This result is quite 
interesting since the nature of the neutrino is controversial, and signals a strong relation 
of this particle with gravity, in an undiscovered way. In the 50's Wheeler |33| had noted 
that the neutrino could be a geometrodynamical particle of some kind (in the next section 
we shall present an Poincare invariance violating extension of the Dirac equation, with 
interesting features nonetheless). 

One outcome of the possible "tachyonicity" of the neutrinos is that, in order to be 
consistent, "tachyonicity" requires the existence of both right handed and left handed 
neutrinos. This is a problem for the neutrino since we have observed only left handed 
neutrinos. Nevertheless, this is already a problem of contemporary neutrino physics, 
which is solved by the seesaw mechanism |34j . 

A negative mass square can be a very strange physical result. However we must note 
that, the mass parameter even in classical mechanics, is a parameter that cannot be 
measured directly even for slow particles [35J. Only energy and momentum are measured 
through interactions and therefore must be real [35]. As noted by the authors of |35j, 
the imaginary mass offends traditional thinking and not observable physics. Nonetheless, 
imaginary mass particles cannot be consistently quantized in the "traditional" quantum 
field theory way, as gravity cannot be quantized too. Perhaps this is another hint that 
makes us speculate about the possible interconnection of gravity and the neutrino field. 

Another quite astonishing feature of tachyonic fermions is that, if such a particle loses 
energy in a medium, it will undergo an acceleration [35J. If we assume that the neutrino 
is a tachyon, there is no direct way for this particle to loose energy, due to the fact that 
it has no electric charge. Nevertheless, it is known that a neutrino travelling through a 
medium can have leading order matter induced magnetic and electric dipole moments [36J . 
The last can be non-zero even for Majorana neutrinos. It is not likely to directly measure 
such effects, however if the neutrino is proven to be tachyon, then these interactions could 
make the neutrino loose energy in various dense medium. This could have observable 
effects coming from astrophysical sources and astrophysical events. Since the energy loss 
of a tachyon causes it's acceleration, then neutrinos loosing energy in a Gamma Ray Burst 
generating background, would eventually accelerate and would reach earth faster than a 
photon. This result is similar to the quantum gravity one, we discussed earlier. 

Relativistic Quantum Mechanics, Localization and A-causality and Superlu- 
minal Propagation 

In this section we shall discuss a curious and rather not so very well known problem that 
appears in the relativistic Dirac equation. We base our analysis in [37]. It is known 
that any solutions of the Dirac equation cannot describe propagation faster than light. 
However in the Foldy-Wouthuysen representation (a unitary transformation of the Dirac 
and other operators, that transforms the Dirac equation into two Klein-Gordon equations), 
superluminal propagation is possible, in terms of a superluminal spreading of the wave 
functions. This problem is very closely related to position operators commuting with the 
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positive energy and in addition to the way we define localization of a particle. If we assume 
a quantum mechanical-consistent definition of the localization of a particle in a subspace 
B of i? 3 , then there is a non-zero probability that this particle can be found outside this 
space B, at a space -Bo. This instantaneous spreading of the wave functions corresponds 
to superluminal propagation. The whole argument can be summarized in terms of the 
following theorem, which states that: 

Let D be a Hilbert space. For every set B <Z B? let F{B) be a bounded self adjoint op- 
erator in D, such that for every function 1(1 belonging in D with (ip\ip) = 1, the probability 
measure {ip\F{B)\ip) satisfies the following properties: 

• {ip\F(B)\ip) = 1 and (<j>\F{B)\<j>) = implies {<f>\if>) = 1 

• There is a self adjoint operator p (the generator of space translations) such that for 
all a e R 3 , 

(e~ l P a Tp\F(B + a)|e-^» = {ip\F(B)\ip) 

Furthermore define a non-constant positive operator H{p). For all non empty sets 
Bi C R? and for all e > 0, there exists a time t e (0,e), for which 



Therefore if F(Bj) represents a position operator, there is a small probability that 
the particle can be outside Bi, to another subset of R? for arbitrarily small time. It is 
expected that in physically reasonable situations, this instantaneous spreading of the wave 
functions must be very small. For a much more detailed analysis we refer the reader to 
reference [37J (page 30). 

An [7(2, 2)-invariant Modification of the Dirac Equation 

The standard Dirac Lagrangian in curved spacetime is built as a gauge theory on a man- 
ifold M, with SL(2, C) structural group. The vielbeins are the sections of the corre- 
sponding principal bundle, the spin connections uj r sll are corrections of the principal bundle, 
and the fields tp r are the sections of the associate bundle with fibre C 4 . A more general 
spinor equation in curved space can be built if we admit the structural group to be the 
total pseudo-unitary group [7(2,2), to which SL(2,C) is an injected subgroup [38]. The 
group C7(2, 2), is used in twistor geometry and conformal field theory. We shall consider a 
very simple form of the modified [7(2, 2)-Dirac Lagrangian (see [38] for details), by taking 
e A = ^A' ^ ne internal metric to be locally Euclidean, that is g^ v = 8^,u- The Lagrangian 
reads, 



(19) 




(20) 



The corresponding equation of motion is, 



Vi^d^ip -Wtp- Ud^d^ip = 



(21) 
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with V,W,U, real constants (some of these have mass dimensions). The above equation 
does not correspond to any relativistic equation of an elementary particle, since no rep- 
resentation of the Poincare group results in equation ()2ip . In addition, the quantization 
of the above equation is not usual in the traditional Quantum Field Theory-way because 
tachyonic solutions appear, and obviously due to the second derivatives in (f2Tj) . Never- 
theless we think that (|2ip describes a particle dynamical system not respecting Poincare 
invariance in 4-dimensions. The reason for considering such an physics- unattractive equa- 
tion, is due to the outcomes, in reference to tachyonic solutions. Indeed, if the solutions 
of (|21|) are of the form, ip = cj)e~ iPl * x , then equation (|2ip . becomes: 

= m0 (22) 

with two solutions for the mass, namely: 

mi, 2 = (^2 ( 2UW + y2± vV 4 + 4UWV 2 ) ) V2 (23) 

Therefore, the solutions to (|21|) . correspond to a superposition of two Dirac waves with 
masses mi^m^- There are values in the parameter space for which the masses can be 
purely imaginary. This means that we may have tachyonic solutions for the wave particle 
■0. In addition there is a mass gap between the two masses. This means that if we choose 
appropriately the parameters, we can obtain a large mass gap between the masses m\ 
and 77i2- If we consider that the field tp could describe the neutrino, both the results we 
described above are quite interesting. Especially the second result, since it is reminiscent 
to the phenomenological outcomes of the seesaw mechanism |34j . Indeed if the two wave 
solutions with masses m\ and mi, describe the left handed and right handed neutrino, 
we can achieve the mass splitting to be huge, so that the right handed neutrino's cannot 
be observed. We must notice that the aforementioned mass splitting occurs even for real 
masses and not only for imaginary masses. This scenario is quite attractive, in regards 
to tachyonic solutions, since tachyons always come in left and right handed components. 
Hence the above mechanism could provide a way to split energetically the two masses, 
keeping one low and the other one decoupled due to it's high mass. 

Conclusions 

In the present paper we studied a toy model of fermions who's main feature is the dy- 
namical generation of a tachyonic mass for the fermions involved. No attempt was made 
to produce phenomenology or imply that the neutrino indeed is a tachyon. The only 
assumption we made is that the toy fermion has a Lagrangian of the form The phe- 
nomenological implications of such a Lagrangian at the Standard Model level are identical 
to those of a massless Dirac fermion. It is quite a mystery why it yields tachyonic mass 
generation when non-perturbative effects are taken into account. 

Furthermore we discussed various implications of a tachyonic neutrino in physical 
processes. We can sum up here the most important ones: 
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• The energy loss of a tachyon results in the acceleration of the particle. Consequently, 
neutrinos loosing energy in a Gamma Ray Burst generating background could reach 
earth faster than a photon. 

• In hot-dense backgrounds, the plasmon decay and the effective neutrino-photon cou- 
pling can be modified and therefore have direct impact on the neutrino emission 
from stars 

• There are no tachyonic neutrinos which are only left handed. 

In addition we provided a toy £7(2, 2) invariant modification of the Dirac equation 
(originally presented in [38j but in a different context) which can admit tachyonic solu- 
tions and can provide a natural seesaw-like mechanism. However this modification we 
provided cannot describe a common sense-acceptable relativistic particle. It is just a 
geometrodynamical toy model with interesting features. 

It is tempting to discuss problems such as the chiral symmetry restoration at finite 
temperature. Nevertheless it is extremely difficult to address these problems without 
knowing how to handle the resulting extrema. The finite temperature effective potential 
of the modified Gross-Neveu model we used is: 



LZ ^^ <v-\-l)\l\ 



k=l 1=0 y 2 



A ((2ir) 2 f-^-\v - iV 



oo 

" E £ "",'„' r T^ ^r-n (25) 

k=l 1=0 



with a = ^jr- and cz\ = ^j^- In addition, 

d + 1 



(26) 



and also, A is a positive integer. In the above, m represents the dynamical generated 
minimum of the effective potential a. This is to be found by minimizing the effective 
potential. Using the zeta regularization technique [55H32] the effective potential can be 
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Z=0 ^ 2 'H" 

We kept the above expression without simplifying in order to have a clear picture of 
the terms appearing. In the case d = 1 poles appear in the expression above, which are 
cancelled using the zeta regularization technique. This can easily be seen Taylor expanding 
around d = 1 + e for e — > 0. 

It is obvious that any question in reference to chiral symmetry restoration, corresponds 
to finding the minimum of the effective potential above, in terms of a = m. The locus 
of local extrema of the thermal effective potential is very complex to handle. Moreover 
we cannot be sure of what an imaginary solution for the thermal mass could indicate, 
therefore in order to avoid un-physical speculations we postpone such an investigation. It 
would be interesting though to study the finite temperature case with the appearance of 
a complex (or not) chemical potential. We hope to address such problems in the future. 
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